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In this paper we investigate the potential of 3D cosmic shear to constrain massive neutrino pa- 
rameters. We find that if the total mass is substantial (near the upper limits from LSS, but setting 
aside the Ly alpha limit for now), then 3D cosmic shear + Planck is very sensitive to neutrino 
mass and one may expect that a next generation photometric redshift survey could constrain the 
number of neutrinos N v and the sum of their masses m v = ^\ mi to an accuracy of AW„ ~ 0.08 
and Am„ ~ 0.03 eV respectively. If in fact the masses are close to zero, then the errors weaken 
to AN V ~ 0.10 and Am„ ~ 0.07 eV. In either case there is a factor 4 improvement over Planck 
alone. We use a Bayesian evidence method to predict joint expected evidence for N v and m v . We 
find that 3D cosmic shear combined with a Planck prior could provide 'substantial' evidence for 
massive neutrinos and be able to distinguish 'decisively' between many competing massive neutrino 
models. This technique should 'decisively' distinguish between models in which there are no massive 
neutrinos and models in which there are massive neutrinos with \N V — 3| >0.35 and m„ > 0.25 eV. 
We introduce the notion of marginalised and conditional evidence when considering evidence for 
individual parameter values within a multi-parameter model. 
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I. INTRODUCTION 

In this paper we will investigate the potential of 3D 
weak lensing to constrain the properties of, and provide 
evidence for, massive neutrinos. The conclusion that neu- 
trinos have mass, and the resolution of the actual masses, 
would have a profound impact on our understanding of 
particle physics and cosmology. 

As a photon travels from a distant galaxy the path 
it takes is diverted, by the presence of large-scale struc- 
ture along the line of sight resulting in the image of any 
galaxy being slightly distorted. This weak lensing effect 
depends on both the details of the matter power spec- 
trum and growth of structure as well as the geometry of 
the observer-lens-source configuration. 3D weak lensing 
combines this weak lensing information with any redshift 
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information available which then allows for evolving ef- 
fects, for example dark energy, to be investigated. Weak 
lensing (see [48[ for a recent review) has been used to 
constrain cosmological parameters including dark energy 
parameters [33[ , measure the growth of structure [f| |43| 
and map the dark matter distribution as a function of 
redshift [6l|. It has been shown that 3D weak lensing 
has the potential to constrain dark energy parameters 
to an unprecedented degree of accuracy using upcoming 
and future surveys e.g. Pan-STARRS [H and DUNE 
In addition to dark energy parameters, galaxy red- 
shift surveys 0, [2(| H2| and weak lensing tomography (in 
which galaxies are binned in redshift) [l|, [l(J, [21( have the 
potential to constrain the total neutrino mass. In this 
paper we consider a novel technique, 3D cosmic shear 
[1, [H, [24], [33| in which galaxies are not binned in red- 
shift and the full 3D shear field is used, thus maximis- 
ing the information extracted. We the report 3D cos- 
mic shear performance in constraining neutrino proper- 
ties and also present Bayesian evidence calculations that 
will show whether future weak lensing surveys will find 
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convincing evidence for massive neutrinos. 

In the standard model of particle physics neutrinos 
must have zero mass by definition; if neutrinos have mass 
this may be a signal of non-standard neutrino interac- 
tions or Higgs mechanisms, evidence that the standard 
model is not renormalisable, or extra (stringy) dimen- 
sions (we refer to neutrino mass and oscillation reviews 
f32l 1371 167| and references therein). In cosmology neutri- 
nos play a role in structure formation by damping struc- 
ture on small scales. They can be categorised as a hot 
component of dark matter, though recent results (for ex- 
ample WMAP [60]) rule out massive neutrinos as the 
dominant dark matter component. Beyond the affect of 
neutrinos on large scale structure they are of further in- 
terest in cosmology since by directly observing neutrinos 
they could provide a window on the early Universe be- 
yond the surface of last scattering to the epoch of elec- 
troweak unification. Theoretically massive (majorana; 
v a = Vu) neutrinos may provide some explanation for the 
baryon asymmetry [l8j . And indeed astronomical consid- 
erations first alluded to neutrino mass since the number 
of neutrinos detected from the Sun was much less than 
the expected number from the Sun's luminosity. 

There is currently a substantial and growing amount of 
evidence that neutrinos have mass (for recent reviews see 
[32I [37l ) . This conclusion has been reached using re- 
sults from large particle physics experiments which have 
found that the oscillation of neutrinos from one flavour (e, 
(i or t) to another is needed to explain the observed data. 
SuperKamiokande [Tjj found that only ~ 1/3 the flux of 
muon neutrinos from cosmic ray collisions in the atmo- 
sphere were observed along the line-of-sight through the 
Earth implying an oscillation of to some other flavour 
with a scale length comparable to the radius of the Earth. 
SNO 3] has observed both the total flux of neutrinos 
from the Sun as well as the flux of electron neutrinos v e 
and found that the 'solar neutrino problem' is resolved 
by postulating that v e oscillate to other flavours in the 
high density environment of the Sun's core by gaining a 
small effective mass. Further evidence for neutrino os- 
cillations comes from nuclear reactors (KamLAND, [TTj ]) 
and neutrino beam experiments (K2K, [3]). 

The observed oscillation of neutrinos is linked to the 
implication that neutrinos have mass via the lepton mix- 
ing matrix U 2 whose elements describe the probability 
for one neutrino flavour to oscillate to another 
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this relates the observed neutrino flavours \u a ) to a hi- 
erarchy of neutrino mass states \vi). In general the ele- 
ments of U can be complex. The standard model has 3 
neutrino flavours v a — v e , v^, v T . The lepton mixing ma- 
trix allows for more than 3 neutrino flavours; if there are 
more than 3 then the extra neutrinos are called "sterile" 
since they would not couple to any standard model elec- 
troweak interaction (i.e. would not be associated with 
electrons, muons or tau). 

Assuming 3 neutrino flavours U can be parameterised 



as a product of three Euler rotations U = R23R13R12 
where each rotation describes how one neutrino mass 
state is coupled to another. The elements of the rotation 
matrices depend on the mass difference Sniij = rrii — rrij 
and an angle 0y on which the probability P(yi — ► Vj) 
depends. Thus neutrino oscillation experiments can only 
measure the relative masses of neutrinos not the absolute 
mass scale. 

Currently world neutrino data are consistent with a 
three-flavour mixing framework (see [lBJ and references 
therein) , parameterised in terms of three neutrino masses 
(TOi,m 2 ,m 3 ) and of three mixing angles (#12, #23, #13); 
plus a possible CP violating phase S. Current con- 
straints are <5m 2 3 ~ 0.05 eV and Smn ~ 0.007 eV. There 
are currently no strong constraints on <5wi3 though up- 
coming experiments, for example T2K 27J should mea- 
sure A6*i3 ~ 0.05. Thus current constraints allow for 
two possible orders of the massive neutrino hierarchy: 
mi < m-i < 7713 or the inverted hierarchy 7713 < mi < m-z. 
There are planned particle physics experiments that will 
measure the absolute mass scale via the beta decay of 
Tritium to constrain m„ e [38l . |45| . l6a ] , for example KA- 
TRIN [H, H§| should reach an accuracy of Arrive ~ 0.35 
eV. There are arguments (see [13]) that the requirement 
on the accuracy of the absolute mass scale needed to 
break the hierarchy or inverted-hierarchy degeneracy is 
Am„ < 0.1 eV. As discussed below this can be achieved 
by cosmological observations. 

Cosmologically massive neutrinos play a role in struc- 
ture formation since free-streaming neutrinos can sup- 
press growth on small scales. Neutrinos streaming from 
an over-density will reduce the amount of matter that 
can gravitationally accumulate by providing an extra ef- 
fective pressure. In the nomenclature of cosmological pa- 
rameter estimation massive neutrinos modify the matter 
power spectrum's growth rate by providing a suppres- 
sion at small scales. It can thus leave key signatures in 
large scale structure data (see, eg. ,UM ) and, to a lesser 
extent, in CMB data (see, e.g.,|41J). Very recently, it 
has also been shown that accurate Lyman-a (Lya) for- 
est data [44|, taken at face value, can improve the current 
CMB+LSS constraints on 77i„ by a factor of ~ 3, with im- 
portant consequences on absolute neutrino mass scenar- 
ios [57j ■ Current cosmological constraints from WMAP 
CMB combined with SDSS BAO and including Lyman- 
a constraints provide an current upper limit on the to- 
tal mass of mj, < 0.42-0.79 eV depending on the as- 
sumptions made, with a median value of m v < 0.67 eV 

[ilESIHilESIsg. 

There has been substantial work in numerically esti- 
mating the growth rate including the presence of massive 
neutrinos (for example [13) HH)- A degeneracy between 
neutrino mass and dark energy parameters arises because 
they both effectively suppress the matter power spectrum 
growth rate as highlighted by [3l| . Optimistically, meth- 
ods which can constrain dark energy parameters well 
should also be able to constrain the neutrino mass pa- 
rameters, and by combining constraints from multiple 
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methods (e.g. CMB, weak lensing) parameter degenera- 
cies should be lifted. 

It is currently believed that the hierarchical mass scale 
of neutrinos implies that the total mass of neutrinos will 
be approximately ~ 0.04 - 0.1 therefore probes which 
are sensitive to this range of values are required to effec- 
tively constrain the neutrino mass. The remainder of 
this paper will highlight the possibility of using 3D cos- 
mic shear to constrain massive neutrino properties. We 
will introduce the methodology and assumptions made 
in Section |TT] and present results in Sections IIIII and IIV1 
in Section [V] we will discuss our conclusions. 



II. METHODOLOGY 

The central quantity in cosmological neutrino mass 
constraints is the fraction of the total matter density that 
is attributed to massive neutrinos /„ = Q, v /Q m which we 
take from [g| (Eqn. 1) 

'-<»» (s&v)Gs&) (2 » 

where m„ = Ylv=o m i * s ^ ne total neutrino mass, a sum 
over all neutrino species each with a mass m^. The ef- 
fect of massive neutrinos on the matter power spectrum 
is commonly expressed using AP(k)/P(k) — [P(k; f u ) — 
P{h;f v = 0)]/P(k;f v — 0) which decreases and is neg- 
ative as the wave-number k increases and power is sup- 
pressed due to the free-streamin g of the neutrinos. We 
use the Eisenstein & Hu (1998) [12[ fitting formula for 
the linear power spectrum which depends on both the 
number of massive neutrino species N v and the total neu- 
trino mass and use the modification of the linear growth 
factor / = d\nS/d\na suggested by [3l| which adds a 
further dependence on /„ and in addition a dependence 
on the dark energy equation of state w = p/p. We use 
the common parameterisation of the dark energy equa- 
tion of state w(a) — wo + w a (l — a) 0. One limitation of 
the Eisenstein & Hu (1998) [12] parameterisation is that 
it assumes that each massive neutrino has the same mass 
i.e. m„ = N v rrii where rrii is the same for all neutrino 
mass eigenstates. Further to this approximation we will 
also treat the number of massive neutrinos as a contin- 
uous parameter that is fitted by data, as opposed to an 
integer number, this can be justified since any light par- 
ticle that does not couple to electrons, photons or ions 
will contribute to the effective number. However these 
assumptions imply that the constraints and predicted ev- 
idence presented are meant to be indicative of the ability 
of 3D cosmic shear to constrain neutrino mass and not as 
entirely representative of the situation as it is currently 
understood - in which there are an integer number of 
mass eigenstates each with a different mass. 

We use a k range of k — 0.001 - 1 Mpc^ 1 for the weak 
lensing Fisher matrix calculations and use the Smith et 
al. (2003) [5l| non-linear correction to the linear power 



spectrum. Note that this is in the quasi-linear regime and 
using wavenumbcrs that are at the limit of the reliability 
of the linear power spectrum fitting formula for massive 
neutrinos, for a recent review of the effect of massive 
neutrinos on the non- linear power spectrum see [54[ . 

A. 3D Cosmic Shear 

The cosmological probes that we will consider in this 
paper are 3D cosmic shear and CMB. We use a CMB 
Planck Fisher matrix which is calculated using CMB- 
FAST (version 4.5.1, [H|) using the method outlined in 
[6l| and will effectively be used as a prior in the results 
presented for 3D cosmic shear. For a general discussion 
of using Planck to constrain neutrino mass parameters 
see [3. 

The weak lensing method we use, 3D cosmic shear 
d, [H, [H, [H| , uses the weak lensing shear and redshift 
information of every galaxy. The 3D shear field is ex- 
panded in spherical harmonics and the covariance of the 
transform coefficients can be used to constrain cosmolog- 
ical parameters. The transform coefficients for a given 
set of azimuthal i and radial k [ftMpc^ 1 ] wave numbers 
are given by summing over all galaxies g; 

fif(Je,£) = ^ 1 a kj l {kr 9 )e- UM9 , (3) 
9 

following the conventions of Q , and we assume a flat sky 
approximation. 6 9 is the angular position of a galaxy on 
the sky, r 9 is the comoving distance to the galaxy and ji 
are spherical Bessel functions. 7 ff is the measured shear of 
the galaxy which parameterises the amount of distortion 
that the galaxy image has obtained due to intervening 
large scale structure. 

This is a novel approach over other 3D weak lensing 
analyses since galaxies are not binned in redshift which 
may cause problems at the bin boundaries and mean in- 
formation loss in averaging over the bins. The binning 
approach, weak lensing tomography, creates a 2D map 
of the galaxies distortions at each redshift and takes the 
cross-correlation between each map to gain some extra 
3D information. Conversely, 3D cosmic shear, presented 
here, uses the entire 3D shear field thus maximising the 
potential for information to be extracted from the galaxy 
image distortions. 

Since the mean of the coefficients in Eqn. ([3]) is zero 
the covariance is varied until it matches that of the data 
[HI, i.e. the covariance is used as the 'signal'. The 3D 
cosmic shear covariance depends on the the lensing ge- 
ometry and the matter power spectrum, so the total pa- 
rameter set that can be constrained is: O m , fide, ^b, h, 
as, wo, w a , n s , the running of the spectral index a n , 
m u and N v . We also include the tensor to scalar ratio 
r = T/S and the optical depth to last scattering r for the 
CMB Fisher matrix calculation. We do not assume spa- 
tial flatness and all results on individual parameters are 
fully marginalised over all other cosmological parameters. 
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Survey 
Area/sqdeg 

^median 

no I sqarcmin 

<7.(z)/(l + z) 



DUNE (fiducial) 
20, 000 
0.90 
35 
0.025 
0.25 



TABLE I: The parameters describing the weak lensing survey 
investigated. 



We use an £ max = 5000 for the 3D cosmic shear analysis 
and a fc max = 1-0 Mpc -1 and use the same assumptions 
presented in [24j . For the Planck constraints we use a 
maximum I of i max — 2500. We will present results for a 
fiducial weak lensing survey which is based on the DUNE 
weak lensing concept, the assumed survey parameters are 
outlined in Table U The z mc dian is the median rcdshift 
of the number density distribution of galaxies with red- 
shift, and no is the observed surface number density of 
galaxies. o~ z (z) describes how the average accuracy with 
which a galaxies position in redshift is known, and er c is 
the statistical variance of the intrinsic observed distor- 
tion of galaxies due to their random orientation. Note 
that we expect the photometric redshift error to have 
a small effect on the predicted statistical constraints as 
shown in [24j], however a good photometric redshift er- 
ror is required to reduce the effect of intrinsic alignment 
systematics as shown in [6| and [34[ . 



B. The Fisher Matrix & Bayesian Evidence 

The results presented in this paper will use the Fisher 
matrix formalism to make predictions of the cosmolog- 
ical parameter errors. The Fisher matrix is defined as 
the second derivative of the likelihood surface about the 
maximum. In the case of Gaussian distributed data with 
zero mean this is given by [H, HiJ [63[ 



F aP = ^Tr [C^CaC^Cp] 



(4) 



where C = S+ N is the theoretical covariance of a partic- 
ular method which consists of signal S and noise N terms. 
The commas in Eqn. (j4|) denote derivatives with respect 
to cosmological parameters about a fiducial cosmology. 
The fiducial cosmology used in this paper is based on the 
WMAP results [gcj|; fi m = 0.27, n de = 0.73, fl b = 0.04, 
h = 0.71, cr 8 = 0.80, w a = -1.0, w a = 0.0, n s = 1.0, 
a n = 0.0. We consider two sets of fiducial value for the 
neutrino mass; one which is in agreement with current 
cosmological constraints, and one in which there are no 
massive neutrinos. This will allow for some discussion 
on the sensitivity of the predicted results to the fiducial 
values. The first m v = 0.66 eV and N v — 3.0 is high 
compared to the expected hierarchical mass scale. We 
justify this since we are using the current constraint on 
the neutrino mass from cosmology. The second set of 
fiducial values are m v = eV and Af„(massless) = 3.0. 



The predicted marginal errors are calculated by tak- 
ing the inverse of the Fisher matrix, the error on the 
i th cosmological parameter is given by A9i = 
To combine constraints from multiple experiments the 
Fisher matrices are summed e.g F to tai = Sensing + Fcub ■ 

In addition to presenting the marginal error on the 
cosmological parameters we will present the expected 
Bayesian evidence that the fiducial survey could achieve 
for massive neutrinos. Computing the evidence allows 
one to distinguish different models rather than constrain 
parameters within a model (e.g. [26j ) as explained be- 
low. A procedure for calculating the expected evidence 
directly from the Fisher matrix was presented in [25j . In 
the case of massive neutrinos there is a natural question 
that an evidence calculation can answer: does the data 
provide evidence for massive neutrinos? Note that this is 
distinct from assuming that massive neutrinos exist and 
using the data to constrain their expected properties. 

The concept of evidence is derived from B ayes' theo- 
rem which relates the probability of model given the data 
p(M\D), to the probability of the data given the model 
p(D\M) 



p(M\D) = 



p(D\M)p(M) 



(5) 



where p(M) is the prior probability on any parameters 
within the model M. We assume two competing models 
M and M'. We also assume that M' is a simpler model 
than M, containing fewer parameters n' < n, and that 
the models are nested i.e. the more complex model M is 
an extension of the simpler model M'. By marginalisa- 
tion p(D\M), known as the evidence, is 



p(D\M)= d9p(D\9,M)p(9\M). 



(6) 



The posterior relative probabilities of the two models, 
regardless of what their model parameters are, is 



p(M'\D) _ p(M') J dff p(D\9',M')p(9'\M') 
p(M\D) ~ p{M) J d9p{D\9,M)p{9\M) ' 



(7) 



By assuming uniform priors on the models, p(M') — 
p(M), this ratio simplifies to the ratio of evidences which 
is called the Bayes Factor, 



Jd9' P (D\9',M')p(9'\M') 
J d9p(D\9,M)p(9\M) ' 



(8) 



It is the evaluation of the Bayes factor that allows one 
to determine whether a data set can distinguish between 
competing models. We wish to forecast whether a future 
survey will be able to distinguish between models. This 
can be done using the Fisher matrix using the following 
expression, given in (25| 



B = (2tt) 



-p/2 V dct F 

Vdet F' 



exp 



1 



1aF a p59[3 



q=l 

(9) 
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with 50 a given by 
56' a = -{F'- l ) a pG K 5ip < a,/3 = l...n,C= l-.-P 

(10) 

where p = n — n' . Note that i* 1 and _F 1 are n x n ma- 
trices, F' is n' x n', and G is an n' x p block of the full 
n x n Fisher matrix F. dtp are the differences in the pa- 
rameters values between models M and M' . AO are any 
prior ranges imposed on the parameters we set AO = 1 
at all times. The expression in Eqn. ([9]), given its im- 
plicit assumption of Gaussian likelihood surfaces, allows 
one to very quickly evaluate the expected evidence. This 
was done in [25j for the case of modified gravity, forecast- 
ing the expected evidence for a single extra parameter 7 
which parameterises any deviation from General Rela- 
tivity. Here we will use Eqn. to calculate the joint 
evidence for two parameters vn v and N v . The Bayesian 
evidence has been extensively studied in cosmology in 
general (for example [H, S, [HE E^l) and in the field of 
dark energy (for example [36l l47l. 153. 158} ) . 

Throughout this paper we will use the Jeffreys [2II ] 
scale in which, InB < 1 is 'inconclusive', 1 < InB < 2.5 
is described as 'substantial' evidence in favour of a model, 
2.5 < \nB < 5 is 'strong', and InB > 5 is 'decisive'. 

III. MARGINAL ERROR RESULTS 

In this Section we will present the predicted marginal 
errors on the massive neutrino parameters could be found 
using 3D cosmic shear in combination with a CMB 
Planck experiment. 

Table [TT] shows the marginalised cosmological parame- 
ter errors predicted for Planck alone and combined with 
the 3D cosmic shear constraints from the fiducial sur- 
vey. We also include the dark energy pivot redshift er- 
ror. The pivot redshift is the point at which the error 
on w p minimises, this is defined by rewriting the dark 
energy parameterisation used w(a) — wq + w a (l — a) to 
w(a) = w p + w a (a p — a) where a p = a p (z p ) can cor- 
respond to any redshift. The Dark Energy Task Force 
(DETF; 0) Figure of Merit (FoM) is defined as being 
proportional to the reciprocal of the area constrained in 
the (w p , w a ) plane at the pivot redshift 

FoM = 1/ Aw p Aw a . (11) 

Note we use the reciprocal of the 1-a two parameter el- 
lipse, the DETF use the 2-a two-parameter ellipse which 
differs from equation (jTT|) by a constant factor. 

By combining 3D cosmic shear constraints with Planck 
the massive neutrino parameters could be constrained 
with marginal errors of Am v ~ 0.03 eV and AN V ~ 0.08 
if the neutrinos are massive, this is a factor of 4 im- 
provement over Planck alone. The dramatic improve- 
ment comes from the lifting of parameter degeneracies 
when the extra constraints are added. As shown in 
[34| without massive neutrinos the fiducial survey de- 
sign could provide a dark energy FoM= 475. The in- 



clusion of the massive neutrino parameters does not de- 
grade this FoM substantially, since the extra parameters 
are well constrained. If there are no massive neutrinos 
then the marginal errors on these parameters degrade, in 
this case the mass and number could be constrained to 
Am„ ~ 0.07 eV and AJV„ ~ 0.10, however this is still a 
factor of 4 improvement over Planck alone. 

This degradation in the marginal error occurs because 
the effect of massive neutrinos on the matter power spec- 
trum and hence on 3D weak lensing is non-linear. If 
the mass of neutrinos is larger then the amount of sup- 
pression at a given scale increases, furthermore the effect 
on the linear growth factor as a function of /„ given by 
[3l| is non-linear. In addition the scale at which power 
is suppressed due to free-streaming varies as a function 
of neutrino mass; neutrinos with lighter masses suppress 
growth at larger scales (higher-k) at a given redshift. [62| 
investigated the effect of the neutrino mass fiducial model 
when making predictions on future marginal errors for a 
galaxy redshift survey and found that if the fiducial value 
°f fv % 0-01 then the marginal error on N v depends 
strongly on the assumed value of /„ (e.g. [12], Figure 
3). One should only expect parameters fiducial values 
to have a small effect on the predicted marginal errors if 
they have a linear affect the covariance of a method i.e. so 
that the derivatives in the Fisher matrix (Eqn. 2]) do not 
change as the point around which the derivative is taken 
changes. The assumed fiducial value of the neutrino pa- 
rameters has a small effect on the errors of most other 
cosmological parameters, given that the fiducial models 
are so different, except for Q m . This is because the con- 
tribution of massive neutrinos to the total mass-energy 
behaves with redshift like a extra matter density. 

Hannestad et al. (2006) ,21] find that, by combin- 
ing the weak lensing tomography constraints from their 
Wide-5 survey (which is similar to our fiducial survey) 
with Planck, Am v = 0.043 eV and AN V = 0.067. They 
assume a massive neutrino fiducial model with m v = 0.07 
eV and N v = 3.0. We find good agreement between this 
result and our massive neutrino fiducial model despite 
using different weak lensing methods, slightly different 
survey designs and different fiducial models. They find 
a Planck-only predicted errors of Am„ = 0.48 eV and 
AN V = 0.19 although they do not include B-modes which 
we do in our CMB Fisher matrix 61], we also checked the 
CMB Fisher matrices with an MCMC analysis and found 
that the error on N v for Planck, assuming that the neu- 
trino mass is fixed, was AN V = 0.084 and for the Fisher 
matrix analysis we find AN V — 0.10 assuming that the 
mass is fixed. 

As shown in [34,] it can be expected that realistically 
the dark energy FoM constraints from 3D cosmic shear 
combined with a Planck prior will be reduced by approx- 
imately a factor of 2 due to photometric, intrinsic align- 
ment and shape measurement systematic effects, this cor- 
respond to a factor of y/2 for each dark energy parameter. 
Since, as highlighted by [3l| . the neutrino mass parame- 
ters affect the power spectrum in a similar way to dark 
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Fiducial values: 


m„ = 0.66 eV & N„ = 


3.00 


m v = eV & N v = 


3.00 


Parameter 


Planck Alone 


3D Cosmic Shear + Planck 


Planck Alone 


3D Cosmic Shear + Planck 


AQ m 


0.0014 


0.0008 


0.0104 


0.0041 


AQdc 


0.0015 


0.0012 


0.0041 


0.0021 


Ah 


0.0167 


0.0055 


0.0148 


0.0060 


Aas 


0.0965 


0.0040 


0.0999 


0.0202 


AQ b 


0.0019 


0.0006 


0.0014 


0.0007 


Awo 


0.5622 


0.0442 


0.6031 


0.0309 


Aw a 


1.8679 


0.2277 


1.9158 


0.1853 


An B 


0.0103 


0.0018 


0.01435 


0.0039 


Aa„ 


0.0083 


0.0044 


0.0074 


0.0046 


Ar 


0.0199 


0.0193 


0.0207 


0.0202 


At 


0.0084 


0.0078 


0.0080 


0.0078 


Am„/eV 


0.2815 


0.0324 


0.3815 


0.0728 


AN V 


0.1144 


0.0836 


0.2807 


0.1042 


Aw p 


0.1177 


0.011 


0.1879 


0.011 


DETF FoM 


5 


400 


3 


490 



TABLE II: Predicted marginalised cosmological parameter errors for Planck alone and combined with the 3D cosmic shear 
constraints from the fiducial survey. We also show the dark energy pivot redshift error and the DETF Figure of Merit, FoM. 
We show results assuming two different sets of fiducial values for the massive neutrino parameters one in which neutrino have 
mass and one in which neutrinos are massless. 



energy one should realistically expect at most a factor 
of a/2 reduction in the combined constraints due to sys- 
tematics. Using this heuristic approximation this yields 
constraints of Am„ ~ 0.04 eV and AN„ ~ 0.11 for a 
massive neutrino fiducial model. 



IV. BAYESIAN EVIDENCE RESULTS 

In the interpretation of the results in the following Sec- 
tion one should keep in mind that the magnitude of the 
Bayes factor shown is a prediction of an experiments abil- 
ity to distinguish one model over another i.e. to what 
level could the experiment in consideration will be able 
to provide evidence for the fiducial model over another 
competing model (or vice versa) where the models are 
distinguished by changes in the parameter values Sifri. 

Note that in the evidence calculation we use the m v — 
eV and N v (massless) =3.0 fiducial values for the Fisher 
matrices since this represents the 'simple' model as de- 
scribed in Section Ul Bl and will allow for statements to be 
made about whether the data could provide evidence for 
a more complicated model containing massive neutrinos 
over a simpler model with no massive neutrinos. 



A. Multi-Parameter Expected Evidence 

Fig. [T] shows the expected evidence contours for m v 
and N u jointly for the fiducial survey design. Note that 
a 5N V — means that N v = 3 and 6m u = means 
that m v = eV i.e. at the fiducial values. The figure 
shows that there is a substantial improvement in combing 
Planck with 3D weak lensing data. On its own Planck 
could only provide at most substantial evidence massive 



neutrinos for models with a large range of massive neu- 
trinos parameter values. The fiducial survey using 3D 
cosmic shear combined with a Planck prior will; 

• Provide substantial evidence for massive neutrinos 
over models in which there are no massive neutri- 
nos, and if the neutrino mass is small 5m u < 0.1 
eV then there will be substantial evidence for these 
models. 

• Be able to decisively distinguish between models in 
which there are no massive neutrinos and models 
in which N v < 3.00 - 0.40 or N„ > 3.00 + 0.40 and 
m v > 0.25 eV. 

• Specifically the experiment could decisively distin- 
guish between models in which there are 3 mass- 
less neutrinos and i) models in which there are 
few N v < 2.6 (possibly zero) massive m u > 0.25 
eV neutrinos ii) models in which there are many 
N v > 3.4 massive m„ > 0.25 eV neutrinos iii) 
models in which there are few N v < 2.6 massless 
m„ = eV neutrinos and iv) models in which there 
are many N v > 3.4 massless m v = eV neutrinos. 

There is a band in which the evidence is inconclusive 
(the black band in Fig. [I}, this represents the boundary 
between where the data would favour the simpler fiducial 
model and the situation in which the data would favour 
a different model (i.e. where the probability of either the 
fiducial or a different model being correct is equal). 

B. Single-Parameter Expected Evidence 

As well as the joint expected evidence on the two mas- 
sive neutrino parameters we can also investigate the ex- 
pected evidence for either parameter individually. When 
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FIG. 1: The expected joint evidence for the number N v and 
mass m„ of neutrinos using 3D cosmic shear and the fiducial 
survey design. White='decisive', lightest gray='substantial', 
darkest gray='strong' and black='inconclusive'. The upper 
panel shows the constraints from Planck alone, the lower 
panel shows the constraints when 3D cosmic shear from the 
fiducial survey is combined with Planck. 

this is done there are two ways in which the other (hid- 
den) parameter(s) can be dealt with; 

• The hidden parameters can either be assumed to be 
fixed at their fiducial values. We will refer to this 
as the conditional evidence. In this case the ex- 
pected evidence presented has the implicit assump- 
tion that the hidden parameters have the value cho- 
sen. 

• The evidence can be integrated over the hidden pa- 
rameters to obtain what we will refer to as the 
marginal evidence. For a multi-parameter model 
which depends on 0i = \,..., n parameters the total ex- 



pected evidence is a function of all these parameters 

B(8i = i „). The marginal evidence on one of the 

parameter 6j is given by 

B(9j) = J ilO , : <l(),. : ..jI0„ Hit), , „:. (12) 

In the case presented here, that of two parameters, 
the total evidence B(N U , m v ) can be integrated to 
obtain the marginal evidence B(N V ) or B(m u ). 

Fig. [J] shows the one-dimensional expected evidence for 
3D cosmic shear combined with a Planck prior. It can 
be seen from both panels in this Figure that by using the 
conditional evidence one can over estimate the evidence 
when the deviation between models is small by upto a fac- 
tor of *~ 5, however when the models being compared are 
very different (large values oiStp) the marginal and condi- 
tional evidences converge. The results drawn from these 
plots are similar to those from the full joint evidence. 3D 
cosmic shear should find substantial evidence that neu- 
trinos are massless if this is the case. Furthermore this 
experiment could decisively distinguish between models 
in which there are no massive neutrinos and models in 
which there are massive neutrinos with m v > 0.25 eV, or 
if the number of neutrinos differs by \N V — 3| > 0.35 

V. CONCLUSION 

In this paper we have shown that 3D cosmic shear has 
the ability to measure the effect that massive neutrinos 
can have on the matter power spectrum, and use this ef- 
fect to place constraints on the total mass of neutrinos 
m v = mi and number of these massive neutrinos N v . 
By combining the results using 3D cosmic shear from 
a next generation photometric redshift survey with the 
constraints from the Planck CMB experiment one could 
expect marginalised errors for the massive neutrino pa- 
rameters of Am„ ~ 0.03 eV and AiV„ ~ 0.08 which is 
a factor of 4 improvement over the constraints using the 
CMB alone. We found that if one assumes in this cal- 
culation that neutrinos are massless then the predicted 
marginal error on these parameters is substantially de- 
graded to Am u ~ 0.07 eV and AiV^ ~ 0.10, however this 
is still an improvement of a factor of 4 over the marginal 
errors from Planck alone using the same fiducial model. 
This increase in the marginal errors occurs because the 
power spectrum is affected by neutrino mass in a non- 
linear way. 

Even by including heuristically including systematic ef- 
fects, using a rule-of-thumb from Kitching et al. (2008), 
the improvement over Planck alone is still a factor of 
3. Comparing with other probes we find that 3D weak 
lensing is competitive; [62] ] find that using a galaxy red- 
shift survey combined with a Planck prior Am„ > 0.025 
eV, [19] find that in combination with Planck constraint 
Lyman-cv experiments could constrain Am u < 0.06 eV. 
However we note that these constraints should be depen- 
dent on the fiducial value of the neutrino mass chosen. 
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We explicitly presented results for a fiducial survey 
which has the characteristics of DUNE however other 
forthcoming surveys are also well suited to do 3D weak 
lensing and should have a similar sensitivity to neutrino 
mass for example Pan-STARRS-1 should yield constraint 
roughly twice that of DUNE [|| Am„ ~ 0.06 eV and 
AN U ~ 0.16. The LSST should yield constraints of 
roughly the same order of magnitude as D UNE. 

Using the expected evidence calculation from [25| we 
have shown that one can expect substantial evidence for 
massive neutrinos if they exist, and furthermore that 
one could decisively distinguish between models in which 
there are no massive neutrinos and models in which there 
are massive neutrinos with \N V — 3| > 0.35 and m v > 0.25 
eV. 

We have introduced the concept of marginal and condi- 
tional evidence and shown that by assuming the value of 
one parameter in a model to be fixed the one-parameter 
evidence can be under or over estimated by upto a fac- 
tor 5. These evidence calculations can be generalised 
to models with an arbitrary number of parameters, and 
the simple application of this algorithm using only the 
Fisher matrix can allow predictions to be made which 
would be prohibitively time consuming using traditional 
evidence calculations (with the caveat that Gaussianity 
is assumed). 

If the constraints predicted in this paper were realised 
then our understanding of massive neutrinos could be 
revolutionised allowing the physics beyond the Standard 
Model which this implies to be understood more entirely. 



FIG. 2: The predicted evidence for the number N v and total 
mass m v of neutrinos individually for 3D cosmic shear using 
the fiducial survey combined with a Planck prior. In each plot 
the solid line show the conditional evidence assuming that 
the other parameter is fixed at its fiducial value, the dashed 
line shows the marginal expected evidence when the possible 
values of the hidden parameter are taken into account, see 
Eqn. (1 1 2 p . The dot-dashed lines show the defining evidence 
limits on the Jeffery's scale where ln£> < 1 is 'inconclusive', 
1 < lnB < 2.5 is 'substantial' evidence in favour of a model, 
2.5 < lnB < 5 is 'strong', and ln_B > 5 is 'decisive'. 
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